AFIT/GEO/ENG /93D-03 


AD-A274  076 

ililltuii 


c 


OIL 

ft 


OPTICAL  WAVELET  TRANSFORM 
FOR 

FINGERPRINT  IDENTIFICATION 
THESIS 

Robert  P.  MacDonald 
Captain,  USAF 

AFIT/GEO/ENG/93D-03 


93-31030 

UIIIIHI 

Approved  for  public  release;  distribution  unlimited 

93  12  22  143 


AFIT/GEO /ENG /93D-03 


OPTICAL  WAVELET  TRANSFORM 
FOR 

FINGERPRINT  IDENTIFICATION 

THESIS 

Presented  to  the  Faculty  of  the  Graduate  School  of  Engineering 
of  the  Air  Force  Institute  of  Technology 
Air  University 
In  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of 
Master  of  Science  in  Electrical  Engineering 


Robert  P.  MacDonald,  B.S.E.E. 

Captain,  USAF 

15  December,  1993 

DWC^ALnTwsPE03BDs 
Approved  for  public  release;  distribution  unlimited 


Preface 


I  can  honestly  say  that  while  working  on  this  thesis,  I  spent  a  good  portion 
of  that  time  in  the  dark,  both  literally  making  holograms  and  figuratively  searching 
for  a  clue.  I  could  not  have  accomplished  anything  without  the  support  of  others.  I 
would  like  to  recognize  the  following  individuals  for  their  guidance.  My  committee 
members,  Dr.  Greg  Warhola  and  Dr.  Dennis  Ruck  for  their  insights.  Dr.  Steve 
Rogers,  who  fortunately  for  AFIT  and  my  research,  decided  to  stay  on  as  a  civilian 
professor  here  at  AFIT.  For  that  matter  I  should  thank  the  Dean  for  making  the 
moves  to  retain  Dr.  Rogers’  talents  here.  Hopefully  many  more  future  students  will 
benefit  from  his  mentoring.  Dr  Byron  Welsh  for  his  enthusiasm  about  my  research. 
And  I  would  especially  like  to  thank  Captain  (soon  to  be  Dr.)  Ken  Fielding  and 
Captain  (all  but  hooded,  Dr.)  Tom  Burns  for  their  combined  wisdom,  without  which 
I  would  have  remained  within  that  caliginous  region  of  unenlightenment. 

This  research  has  been  a  great  experience,  one  that  I  will  always  remember 
fondly.  There  were  various  moments  of  despair,  but  with  the  support  from  the 
“back  of  the  lab  club” ,  specifically  Gary  Shartle,  Kim  McCrae,  John  Keller,  Neale 
Prescott,  Curtis  Martin,  and  Martin  Chin,  those  moments  were  bearable.  I  would 
also  like  to  thank  Jo  Jo  for  her  continued  encouragement  and  wish  her  good  fortune 
in  her  scholastic  endeavors. 

MAC 


Robert  P.  MacDonald 


Table  of  Contents 


Page 

Preface  .  ii 

List  of  Figures  .  v 

Abstract .  viii 

I.  Introduction .  1 

1.1  Summary  of  Current  Knowledge .  1 

1.2  Problem  Statement .  3 

1.3  Scope  .  3 

1.4  Outline  of  Thesis .  3 

II.  Literature  Review .  5 

2.1  Background .  5 

2.2  Transformation  Techniques .  5 

2.2.1  Fourier  Methods  .  5 

2.2.2  Wavelet  Methods .  6 

2.2.3  Analytical  Comparison . 8 

2.2.4  Experimental  Comparison .  8 

2.2.5  Symmetric  Wavelet  Transform .  9 

III.  Methodology  .  12 

3.1  Wavelet  Filter .  12 

3.1.1  Mathematical  Development  of  Wavelet  and  Basis 

Filters .  13 

3.1.2  Binary  Detour  Phase  Hologram .  16 

3.1.3  Two  Dimensional  Filters .  20 


iii 


Page 

3.2  Optical  Architecture .  26 

3.3  Optical  Wavelet  Fingerprint  Identification .  29 

3.3.1  Vander  Lugt  Filter .  29 

IV.  Results .  32 

4.1  Characterization  of  Optical  Biorthogonal  Wavelet  TVansform  32 

4.1.1  Input  Image  of  Letter  “E” .  32 

4.1.2  Horizontal  Filter  Results .  32 

4.1.3  Vertical  Filter  Results .  36 

4.1.4  Diagonal  Filter  Results .  36 

4.2  Optical  Biorthogonal  Wavelet  Transform  of  Lenna .  41 

4.3  Optical  Biorthogonal  Wavelet  lYansform  of  a  Fingerprint  .  42 

4.4  Improved  Optical  Wavelet  Fingerprint  Identifier .  51 

4.5  Test  for  Rotational  Invariance .  53 

4.6  Summary .  59 

V.  Conclusions .  60 

5.1  Future  Applications .  60 

5.2  Recommendations .  61 

Appendix  A.  C  Program  used  to  Create  Detour  Phase  Computer  Gener¬ 
ated  Hologram .  62 

Bibliography .  66 

Vita .  68 


List  of  Figures 


Figure  Page 

1.  Two  Dimensional  Decomposition  Algorithm  (1) .  7 

2.  One  Dimensional  Two-Channel  Subband  Coder .  10 

3.  Diagram  of  Optical  Correlator  .  12 

4.  Basis  function  <j>(x)  and  its  Fourier  spectrum  <f>(x) .  15 

5.  Wavelet  function  ip(x)  and  the  magnitude  of  its  Fourier  spectrum  ij>(x)  15 

6.  Sample  Cell  of  Detour  Phase  Hologram  (30) .  17 

7.  Transparency  of  Detour  Phase  Hologram .  21 

8.  Dekacon  Camera  System .  22 

9.  Horizontal  Wavelet  Filter:  <t>(x)rf;(y) .  23 

10.  Vertical  Wavelet  Filter:  ip(x)<f>(y) .  24 

11.  Diagonal  Wavelet  Filter:  ^(x)^(y) .  25 

12.  Adjustable  Scale  Optical  Correlator .  26 

13.  Improved  Optical  Wavelet  Fingerprint  Identification  Arrangement  .  .  28 

14.  Vander  Lugt  Filter  (18) .  29 

15.  HC-500  Thermoplastic  Recording  Process  (24) .  31 

16.  Capital  Letter  UE“ .  33 

17.  Optical  Wavelet  Transform  for  Horizontal  Filter  on  Letter  E,  m  =  0 

Resolution  .  34 

18.  Digital  Wavelet  Transform  for  Horizontal  Filter  on  Letter  E,  m  =  0 

Resolution  .  34 

19.  Optical  Wavelet  Transform  for  Horizontal  Filter  on  Letter  E,  m  =  1 

Resolution  .  35 

20.  Digital  Wavelet  Transform  for  Horizontal  Filter  on  Letter  E,  m  =  1 

Resolution  .  35 

21.  Optical  Wavelet  Transform  for  Vertical  Filter  on  Letter  E,  m  —  0  Res¬ 
olution  .  37 


v 


Figure 

22. 

23. 

24. 

25. 

26. 

27. 

28. 

29. 

30. 


31. 


32. 

33. 

34. 

35. 

36. 

37. 


Page 


Digital  Wavelet  TVansform  for  Vertical  Filter  on  Letter  E,  m  =  0  Res¬ 
olution  .  37 

Optical  Wavelet  TVansform  for  Vertical  Filter  on  Letter  E,  m  =  1  Res¬ 
olution  .  38 

Digital  Wavelet  TVansform  for  Vertical  Filter  on  Letter  E,  m  =  1  Res¬ 
olution  .  38 


Optical  Wavelet  TVansform  for  Diagonal  Filter  on  Letter  E,  m  =  0 
Resolution  .  39 

Digital  Wavelet  TVansform  for  Diagonal  Filter  on  Letter  E,  m  =  0 
Resolution  .  39 


Optical  Wavelet  TVansform  for  Diagonal  Filter  on  Letter  E,  m  =  1 
Resolution  .  40 

Digital  Wavelet  TVansform  for  Diagonal  Filter  on  Letter  E,  m  =  1 
Resolution  .  40 

Lenna .  41 


Horizontal  Wavelet  TVansform  Results  for  Lenna:  (a)  m  =  0  (b)  m  =  1 
(c)  m  =  2 .  43 

Vertical  Wavelet  Transform  Results  for  Lenna:  (a)  m  =  0  (b)  m  =  1 
(c)  m  =  2 .  44 

Diagonal  Wavelet  TVansform  Results  for  Lenna:  (a)  m  =  0  (b)  m  =  1 
(c)  m  =  2 .  45 

Fingerprint .  46 

Horizontal  Wavelet  Transform  Results  for  Fingerprint:  (a)  m  =  0  (b) 
m  =  1  (c)  m  =  2 .  48 

Vertical  Wavelet  TVansform  Results  for  Fingerprint:  (a)  m  =  0  (b) 
m  =  1  (c)  m  =  2 .  49 

Diagonal  Wavelet  TVansform  Results  for  Fingerprint:  (a)  m  =  0  (b) 
m  =  1  (c)  m  =  2 .  50 

Improved  Optical  Wavelet  Fingerprint  Identifier  Results:  Discrimina¬ 
tion  Against  Subject  1  MSF .  52 


vi 


Figure  Page 


38.  Auto-Correlatiou  Peak  vs  Offset  Rotation  for  Normal  and  Wavelet 

Transformed  Fingerprint  Image  at  Various  Rotations .  54 

39.  Comparison  of  (a)  Five  Degree  Offset  Wavelet  Auto-correlation  and  (b) 
Non-Rotated  Subject  Two  Wavelet  Cross-correlation  with  Subject  One  55 

40.  MSF  Correlation  Peak  Plot  for  Normal  Non-Rotated  Fingerprint  .  .  56 

41 .  MSF  Correlation  Peak  Plot  for  Wavelet  Transformed  Non-Rotated  Fin¬ 
gerprint  .  56 

42.  MSF  Correlation  Peak  Plot  for  Normal  3  Degree  Offset  Fingerprint  .  57 

43.  MSF  Correlation  Peak  Plot  for  Wavelet  Transformed  3  Degree  Offset 

Fingerprint .  57 

44.  MSF  Correlation  Peak  Plot  for  Normal  5  Degree  Offset  Fingerprint  .  58 

45.  MSF  Correlation  Peak  Plot  for  Wavelet  Transformed  5  Degree  Offset 

Fingerprint .  58 


vn 


AFIT/GEO/ENG/93D-02 


Abstract 

The  Federal  Bureau  of  Investigation  (FBI)  has  recently  sanctioned  a  wavelet 
fingerprint  image  compression  algorithm  developed  for  reducing  storage  requirements 
of  digitized  fingerprints.  This  research  implements  an  optical  wavelet  transform  of 
a  fingerprint  image,  as  the  first  step  in  an  optical  fingerprint  identification  pro¬ 
cess.  Wavelet  filters  are  created  from  computer  generated  holograms  of  biorthogonal 
wavelets,  the  same  wavelets  implemented  in  the  FBI  algorithm.  Using  a  detour  phase 
holographic  technique,  a  complex  binary  filter  mask  is  created  with  both  symmetry 
and  linear  phase.  The  wavelet  transform  is  implemented  with  continuous  shift  using 
an  optical  correlation  between  binarized  fingerprints  written  on  a  Magneto-Optic 
Spatial  Light  Modulator  (MOSLM)  and  the  biorthogonal  wavelet  filters.  A  tele¬ 
scopic  lens  combination  scales  the  transformed  fingerprint  onto  the  filters,  provid¬ 
ing  a  means  of  adjusting  the  biorthogonal  wavelet  filter  dilation  continuously.  The 
wavelet  transformed  fingerprint  is  then  applied  to  an  optical  fingerprint  identification 
process.  Comparison  between  normal  fingerprints  and  wavelet  transformed  finger¬ 
prints  shows  improvement  in  the  optical  identification  process,  in  terms  of  rotational 
invariance. 
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OPTICAL  WAVELET  TRANSFORM 
FOR 

FINGERPRINT  IDENTIFICATION 


I.  Introduction 

Automatic,  real-time,  personal  identification  is  an  important  capability  in  a 
security  conscious  environment.  In  an  age  of  almost  unlimited  access  of  informa¬ 
tion  through  computer  networks,  controlled  access  through  automated  recognition 
could  maintain  control  with  limited  oversight  of  facilities.  Whether  for  securing 
access  to  vital  areas  or  assuring  proper  use  of  computer  systems,  this  technology 
is  a  requirement.  Law  enforcement  agencies  use  the  uniqueness  of  the  fingerprint 
for  identifying  and  tracking  criminals.  The  fingerprint  remains  as  one  of  the  con¬ 
stants  against  which  personal  identity  can  be  matched.  A  security  system  can  take 
advantage  of  this  uniqueness  for  controlling  access.  Optics  provide  a  method  of  im¬ 
plementing  a  fingerprint  identification  scheme  in  real-time.  This  research  affects  an 
optical  fingerprint  identifier  with  improvement  in  terms  of  rotational  invariance. 

1.1  Summary  of  Current  Knowledge 

An  optical  correlation  is  an  effective  method  for  discriminating  fingerprints. 
This  procedure  has  been  implemented  real-time  using  a  90-degree  prism  and  a 
Magneto-Optic  Spatial  Light  Modulator  (MOSLM)  (12).  The  advantage  of  this 
configuration  is  the  compactness  of  the  design  stemming  from  a  1-/ joint  transform 
correlator.  The  MOSLM  presents  both  the  input  fingerprint  read  in  on  the  prism 
and  the  print  to  be  matched  against,  making  the  design  easy  to  implement  in  a 
compact  portable  form.  However,  this  configuration  is  very  sensitive  to  rotational 
alignment  of  the  input  print  on  the  prism.  It  has  been  shown  that  the  correlation 
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peak  for  fingerprint  images  drops  off  greater  than  3  dB  for  alignment  offsets  greater 
than  ±1.5°  (14).  This  is  due  to  the  limited  rotational  symmetry  in  the  fingerprint. 
Extraction  of  rotationally  symmetric  fingerprint  features  would  help  to  alleviate  this 
problem,  as  long  as  those  features  still  contain  enough  uniqueness  for  discrimination. 

The  uniqueness  of  the  fingerprint  image  lies  in  the  construction  of  the  ridges. 
Law  enforcement  agencies  single  out  points  where  the  ridges  intersect  and  label  them 
minutiae  points.  These  minutiae  points  are  then  used  for  identifying  the  owner  of 
the  fingerprint.  However,  the  uniqueness  of  the  fingerprint  is  not  limited  to  the  ridge 
layout.  There  are  finer  details  which  are  as  distinctive  as  the  minutiae  points;  the  size 
and  shape  of  the  pores  for  instance  (16).  Since  pore  shapes  are  circular,  they  can  be 
considered  rotationally  symmetric.  An  image  processing  routine  which  emphasized 
those  symmetric,  unique  features  would  contribute  towards  circular  symmetry. 

One  such  routine  is  the  wavelet  transform.  The  wavelet  transform  translates 
signals  into  a  time(space)-frequency  representation.  Representing  image  data  using 
the  wavelet  transform  allows  localization  of  image  spectrum  content.  The  wavelet 
transform  is  readily  implemented  through  correlation  of  the  input  image  and  a 
wavelet  filter.  For  this  reason,  the  wavelet  transform  has  become  a  useful  tool  in  op¬ 
tical  image  processing  (27, 13, 21) .  The  wavelet  transformation  of  a  two-dimensional 
input  function,  f(x,  y)  is  represented  by: 

Wkf(a,b,C,d)=  J  Jeffry)  (1) 

Where  h(£^£,  kf^)  represents  the  two  dimensional  scaled  and  shifted  “mother” 
wavelet. 

The  wavelet  transform  has  been  used  for  detection  of  objects  in  a  2-D  scene 
(10).  Because  the  wavelet  transform  is  an  effective  edge  detector,  it  suppressed 
background  clutter  and  emphasized  the  target  image’s  edges.  The  fingerprint  image 
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is  merely  a  series  of  “edges” .  In  the  same  manner  the  wavelet  transform  can  be  used 
for  emphasizing  the  rotationally  symmetric  details  of  a  fingerprint  image. 

1.2  Problem  Statement 

This  thesis  implements  an  optical  wavelet  transform  on  a  fingerprint  image  to 
enhance  the  identification  process. 

1.3  Scope 

This  research  augments  fingerprint  identification,  accomplishable  in  two  phases. 
The  first  phase  optically  computes  the  wavelet  transform  of  images  for  storage.  Re¬ 
sults  of  the  optical  wavelet  transform  are  compared  to  results  of  previous  research. 
The  second  phase  is  an  optical  correlation  between  wavelet  transformed  fingerprint 
images  for  identification.  The  wavelet  transformed  fingerprint  is  examined  for  iden¬ 
tification  uniqueness  and  rotational  tolerance.  This  research  does  not  explore  al¬ 
ternative  methods  of  matching  fingerprints.  There  is  no  attempt  to  perform  feature 
extraction  on  the  fingerprint  image  for  classification  through  standard  pattern  recog¬ 
nition  algorithms.  Feature  extraction,  a  process  of  identifying  those  features  which 
make  the  image  unique,  emphasizes  a  different  research  area. 

1-4  Outline  of  Thesis 

This  document  is  organized  as  follows.  Chapter  II  provides  the  background 
which  lead  to  the  decision  to  focus  this  research  on  fingerprint  analysis.  The  sym¬ 
metric  biorthogonal  wavelet  transform  is  discussed  as  the  objective  of  this  research. 
Chapter  III  presents  the  various  technologies  involved  in  this  research.  Two  methods 
of  creating  a  hologram  are  presented  in  detail.  The  results  are  presented  in  Chapter 
IV.  The  optical  biorthogonal  wavelet  transform  is  characterized  using  three  images 
for  transformation.  Results  from  an  optical  identifier  are  also  presented.  Finally, 
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Chapter  V  presents  conclusions  and  recommendations  for  future  research  in  this 
area. 
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II.  Literature  Review 


2.1  Background 

Identification  through  fingerprints  is  a  major  tool  in  law  enforcement.  In  order 
to  make  a  positive  identification  of  a  suspect,  that  subject’s  fingerprint  must  be 
on-file  for  comparison  to  the  imprints  left  at  a  crime  scene.  The  digital  storage 
of  each  individual  fingerprint  requires  about  9.8  million  bytes  of  data,  resulting  in 
approximately  245  trillion  bytes  for  storage  of  all  the  fingerprint  cards  on-file  at 
the  FBI  (15).  This  number  can  be  lowered  by  compressing  the  fingerprint  image  in 
preparation  for  storage. 

Image  compression,  which  reduces  the  data  storage  requirement,  is  a  three- 
stage  process:  invertible  transformation;  quantization;  and  redundancy  removal  (31). 
The  transformation  stage,  where  the  image  coefficients  are  translated  to  a  new  co¬ 
ordinate  system,  is  the  most  important  component  of  the  algorithm  (31).  Although 
the  transformation  algorithm  requires  complex  mathematics  to  implement,  its  ben¬ 
efits  are  numerous,  especially  when  the  algorithm  is  accomplished  in  real-time.  This 
literature  review  discusses  four  transformation  techniques  and  a  novel  real-time  im¬ 
plementation. 

2.2  Transformation  Techniques 

2.2.1  Fourier  Methods.  Traditional  methods  for  transformation  are  based 
on  the  Fourier  transformation  process.  The  Fourier  kernel,  equation  2,  an  exponen¬ 
tial,  is  the  basis  set  for  the  traditional  methods. 

e-i2ir(/,*+/r»)  (2) 

Two  techniques  which  are  used  in  image  processing  are  distinguished  by  the  segmen¬ 
tation  of  the  image  and  processing  on  that  segment.  The  standard  transformation 
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technique  for  picture  compression,  established  by  the  Joint  Photographic  Experts 
Group  (JPEG),  is  a  discrete  cosine  transform  on  8  X  8  pixel  blocks  (31).  JPEG 
is  an  effective  means  of  compressing  image  data  and  is  considered  the  “application 
to  beat  in  most  applications”  (31).  An  alternative  to  JPEG  is  the  Localized  Co¬ 
sine  TVansform  (LOT).  Developed  in  an  attempt  to  solve  the  blockiness  problem  of 
JPEG,  the  LCT  technique  has  overlapping  segments  of  bell-shaped  functions  instead 
of  the  discrete  rectangle-shaped  functions  associated  with  JPEG  (31). 

2.2.2  Wavelet  Methods.  A  more  recent  development,  the  wavelet  transfor¬ 
mation  uses  a  different  basis  set  for  its  transformation,  the  wavelet  basis.  Several 
wavelet  functions  have  been  developed  (29).  Because  there  are  several  basis  functions 
for  the  wavelet  transformation,  one  key  advantage  of  this  method  over  traditional 
Fourier  methods  is  the  flexibility  in  the  basis  selection.  The  Orthogonal  Wavelet 
Basis  (OWB)  is  distinguishable  from  the  Best-Basis  Wavelet  (BBW)  in  that  the 
Best-Basis  algorithm  bases  its  selection  for  the  basis  function  on  performance.  The 
BBW  compares  several  transformation  outputs.  The  OWB  is  a  simpler  algorithm 
which  uses  only  one  basis  for  its  algorithm  (31). 
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Figure  1.  Two  Dimensional  Decomposition  Algorithm  (1). 


For  simplicity  and  speed,  a  multiresolution  analysis  (MRA)  scheme  was  de¬ 
veloped  by  Mallat  for  computing  the  wavelet  transform  coefficients  at  dyadic  values 
of  the  scale  and  shift  variables  (23).  The  MRA  is  a  sub-band  coding  scheme  us¬ 
ing  a  quadrature  mirror  filter  bank  for  analysis  and  perfect  reconstruction,  where 
each  level  of  decomposition  filters  the  input  image  with  a  dyadically  scaled  wavelet. 
Figure  1  represents  one  stage  of  the  multiscale  decomposition  algorithm,  where  the 
wavelet  coefficients  are  captured  in  the  horizontal,  vertical,  and  diagonal  detail  im¬ 
ages  at  each  resolution  level  (1).  In  the  flow  diagram  in  figure  1,  the  input  image  is 
fed  to  two  separate  filters,  hO  and  hi,  which  are  effectively  low  pass  and  band(high) 
pass  filters  respectively.  Notice  that  in  the  decomposition  algorithm  each  channel 
is  downsampled  by  the  dyadic  value  which  maintains  the  proper  shift  and  scaling 
operation  for  the  next  level  of  decomposition.  The  filter  operation  is  accomplished 
in  both  dimensions,  rows  and  columns,  creating  four  filtered  versions  of  the  original 
input.  The  three  channels  labeled  details  in  figure  1  are  stored  as  the  wavelet  trans¬ 
form  coefficients,  with  the  fourth  channel  returned  to  the  input  stage  for  the  next 
level  of  decomposition. 


Each  level  of  decomposition  scales  the  filters,  resulting  in  coarser  approximation 
and  detail  coefficients.  In  equations  3  and  4,  the  scaling  and  wavelet  functions,  <f> 
and  ip,  are  scaled  by  the  variable  m  which  corresponds  to  a  level  of  decomposition. 
The  first  level  of  decomposition  is  considered  the  m  =  0  level  of  decomposition,  with 
each  succesive  level  represented  by  m  =  1,2,  •  •  •. 


2~m/2<p(2~mt  -  n) 

(3) 

2~m/2ip(2~mt  -  n) 

(4) 

2.2.3  Analytical  Comparison.  The  wavelet  transformation  has  advantages 
over  the  traditional  Fourier  based  techniques.  It  is  capable  of  representing  the  image 
at  several  scales  in  a  redundant  manner,  due  to  the  multi-scale  property  of  the 
wavelet  transformation  (27).  This  capability  is  useful  for  extracting  characteristics  of 
a  fingerprint  image.  The  wavelet  transform  displays  larger  scale  characteristics,  such 
as  the  fingerprint  ridges,  as  well  as  smaller  scale  characteristics,  such  as  the  position 
and  shape  of  the  pores.  It  is  important  to  access  both  levels  of  detail  (16).  Since 
prints  from  a  crime  scene  are  usually  less  than  perfect,  investigators  must  extract 
as  much  information  as  possible  from  the  available  print  to  match  with  existing 
prints  from  a  data  base.  The  wavelet  transform  can  also  specify  exact  location  of 
frequency  spectrum.  The  localization  of  the  image  in  both  time  and  frequency  makes 
the  wavelet  transformation  more  efficient  for  analysis  of  digital  images  (29). 

2.2.4  Experimental  Comparison.  The  FBI  directed  a  study  which  com¬ 
pared  four  transformation  methods  for  compressing  fingerprint  images:  JPEG;  LCT; 
OWB;  and  BBW  (31, 15).  According  to  the  FBI,  there  are  two  potential  users  of  the 
compressed  fingerprint  image:  automatic  recognition  algorithms  and  human  experts. 
These  users’  analyses  of  the  compressed  image  established  the  comparison  basis  for 
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the  transformation  methods.  When  read  by  the  recognition  algorithms,  the  com¬ 
pressed  image  is  expected  to  contain  the  same  minutiae  relationships  as  the  original 
image.  Additionally,  the  compressed  image  should  stand  up  to  expert  review  for 
detail  (15,  16). 

The  wavelet  transformation  outperformed  the  Fourier  methods  in  both  areas 
of  comparison.  The  compressed  images  were  fed  into  the  Home  Office  Matcher,  an 
algorithm  employed  by  the  FBI  for  matching  suspect  fingerprint  images  to  stored 
ones.  Images  compressed  by  the  wavelet  transformation  resulted  in  an  increased 
number  of  correctly  detected  aspects  compared  to  the  original  images,  while  the 
images  compressed  by  the  JPEG  and  LCT  transformations  resulted  in  a  reduction 
in  detected  aspects  (15).  The  wavelet  compressed  images  also  stood  up  better  to 
close  inspection  by  fingerprint  experts.  When  the  compressed  images  were  brought 
under  close  inspection,  the  images  compressed  by  the  wavelet  transform  showed 
minimal  distortion  at  20:1  compression  ratio,  while  images  compressed  by  the  Fourier 
techniques  showed  excessive  distortion  at  that  same  compression  ratio  (16).  Overall, 
the  wavelet  transformations  performed  considerably  better  than  the  Fourier-based 
techniques  when  the  compression  ratio  exceeded  10:1  (31,  16). 

2.2.5  Symmetric  Wavelet  Transform.  The  standard  for  fingerprint  com¬ 
pression  adopted  by  the  FBI  involves  a  symmetric  wavelet  transformation  (4).  The 
symmetric  wavelet  transform  takes  advantage  of  symmetric  extension  of  the  finite 
input  signal  (fingerprint  image),  accounting  for  the  boundary  conditions  in  a  manner 
which  eliminates  •  ae  jump  discontinuity  introduced  by  a  periodic  extension  (5).  A 
symmetric  extension  has  advantage  over  periodic  extension  in  terms  of  compressibil¬ 
ity. 

In  image  coding,  quantization  errors  will  often  be  most  prominent  around 
the  edges  in  the  images;  it  is  a  property  of  our  visual  system  that  we 
are  more  tolerant  of  symmetric  errors  than  asymmetric  ones.  In  other 
words,  less  asymmetry  would  result  in  greater  :ompressibility  for  the 
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same  perceptual  error.  Moreover,  symmetric  filters  make  it  easier  to  deal 
with  the  boundaries  of  the  image  (17). 


Input 

Image 


Figure  2.  One  Dimensional  Two-Channel  Subband  Coder 


The  decomposition  of  a  signal  into  wavelet  coefficients  can  be  compared  to  a 
Quadrature  Mirror  Filter  (QMF)  in  structure.  In  order  to  use  symmetrical  filters 
and  maintain  linear  phase,  it  is  necessary  that  the  analysis  filters  differ  from  the 
synthesis  filters  (9).  An  input  image  is  decomposed  using  the  filters  hO  and  hi  (see 
figure  2),  fed  to  a  signal  processing  algorithm  (i.e.  image  compression  through  vector 
quantization),  and  recreated  using  a  different  set  of  filters,  fD  and  fl.  The  coefficients 
for  the  synthesis  filters  are  generated  from  the  analysis  filter  coefficients.  The  filter 
coefficients  for  fO  and  fl  are  related  to  the  analysis  filters  hO  and  hi  by  the  following 
relationships: 


hl„  = 


fl.  =  (-1)"+1M)_„+, 


Id  this  manner,  the  filter  bank  maintains  linear  phase  and  symmetry. 
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Thus  the  wavelet  transform  is  more  effective  at  compressing  fingerprint  data 
than  traditional  Fourier  methods,  displaying  improvement  in  both  distortion  analy¬ 
sis  and  recognition.  The  same  symmetric  biorthogonal  wavelet  transform  could  be 
applied  to  the  fingerprint  for  image  analysis,  enhancing  the  identification  process. 
The  technologies  involved  in  implementing  the  optical  biorthogonal  wavelet  trans¬ 
form  are  discussed  in  the  next  chapter.  The  wavelet  transform  is  then  applied  to  an 
optical  fingerprint  identifier. 
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III.  Methodology 


This  research  combines  two  holographic  techniques  for  optical  fingerprint  iden¬ 
tification.  This  Chapter  describes  the  methods  used  to  construct  and  implement  the 
holograms.  First,  the  biorthogonal  wavelet  filter  used  by  the  FBI  for  image  compres¬ 
sion  is  constructed.  The  wavelet  filter  is  the  Fourier  spectrum  representation  of  the 
wavelet  function.  The  wavelet  filter  was  created  using  a  binary,  detour  phase  tech¬ 
nique.  The  wavelet  transform  was  realized  with  continuous  shift  variables  through 
an  optical  correlation  between  the  input  image  and  the  wavelet  filter.  A  telescopic 
lens  arrangement  adjusted  the  scale  variable  of  the  wavelet  transform  continuously. 
The  second  hologram  was  used  as  a  matched  spatial  filter  (MSF)  created  using  the 
Vander  Lugt  method  with  a  reference  beam  and  a  thermoplastic  film  plate. 

3.1  Wavelet  Filter 

The  wavelet  transform  is  simply  a  correlation  between  the  input  image  and 
the  wavelet  function  (Equation  1).  Optics  implements  the  wavelet  transform  in 
real-time  with  continuous  shift  variables  using  optical  correlation.  A  collimated 
laser  illuminates  the  input  image,  the  beam  is  passed  through  a  positive  lens,  a 
representation  of  the  wavelet  filter  is  positioned  in  the  back  focal  plane  of  the  positive 
lens,  and  the  resulting  image  is  passed  through  a  second  lens  resulting  in  the  wavelet 
transformed  image.  This  arrangement  is  shown  in  Figure  3. 

Wavelet  Wavelet 

Input  niter  Transform 

Plane  LI  Plane  LI  Plane 


Figure  3.  Diagram  of  Optical  Correlator 
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The  wavelet  filter  was  a  computer  generated  hologram,  representing  the  Fourier 
transform  of  the  wavelet  function.  This  complex  mathematical  representation  was 
constructed  from  the  finite  impulse  response  (FIR)  coefficients  for  the  wavelet. 

3.1.1  Mathematical  Development  of  Wavelet  and  Basis  Filters.  This  re¬ 
search  implemented  the  same  wavelet  function  chosen  by  the  FBI  for  image  com¬ 
pression.  The  symmetric  FIR  filters  used  in  the  creation  of  the  biorthogonal  wavelet 
filter  were  chosen  because  of  their  demonstrated  ability  to  compress  FBI  fingerprint 
images  (31,  15).  The  values  for  the  nine/seven  tap  symmetric  filter  coefficients  used 
in  this  research  are  given  by  Table  1 . 


Tap 

Exact  Value 

Approximate  Value 

h0o 

-by/2  xj(48  |  X2  |2  -16flx2  +  3)/32 

+0.85269867900940 

hO±i 

-5y/2  *i  (8  |x2  |2  -Rx2  +  3)/8 

+0.37740285561265 

h0±2 

-5>/2  xi(4  |  x2  |2  +4/?x2  -  1)/16 

-0.11062440441842 

h0±3 

— 5\/2  xi(f?x2)/8 

-0.023849465019380 

h0±4 

— 5\/2  xi/64 

+0.037828455506995 

hi -i 

y/2  (6xj  —  l)/16xi 

+0.78848561640566 

hi -2,0 

— y/2  (16xi  —  l)/64xi 

-0.41809227322221 

hl-3,1 

y/2  (2xj  +  l)/32xi 

-0.040689417609558 

hi -4,2 

— \/2  /64xi 

+0.064538882628938 

Table  1.  Analysis  Wavelet  Filter  Coefficients  (5). 
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where  (6): 


x 2 

A 

B 


The  equations  5-8  demonstrate  how  to  generate  the  frequency  domain  representa¬ 
tion  of  each  of  the  filters.  First,  the  transfer  functions  for  each  of  the  filters  were 
generated  from  their  impulse  response. 


">»(*)  = 


"•>(*)  =  ^  (6) 

The  basis  function  in  the  frequency  domain  for  the  analysis  filter  was  determined  by 
equation  7. 

4>{x)  =  -4=  II  mo(2"i*)  (7) 

v2tt  J=1 

The  transform  of  the  wavelet  function  for  the  analysis  filter  was  determined  by  the 
product  of  the  transfer  function  for  the  synthesis  filter  and  the  transformed  basis 
function  evaluated  at  half  the  frequency  as  shown  in  equation  8. 
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ip(x)  = 


(8) 


The  filter  ho  in  Figures  1  and  2  was  represented  by  the  basis  function  in 
equation  7.  The  time  and  frequency  domain  representation  for  the  basis  function 
are  shown  in  Figure  4.  Notice  the  low  pass  characteristic  of  the  basis  filter  0(f)  which 


Figure  4.  Basis  function  4>{x)  and  its  Fourier  spectrum  0(f) 


is  both  real  and  symmetric.  The  filter  hi  in  Figures  1  and  2  was  represented  by  the 
wavelet  function  in  equation  8.  The  time  and  frequency  domain  representation  for 
the  wavelet  function  are  shown  in  Figure  5.  The  wavelet  filter  0(f)  was  a  bandpass 


Figure  5.  Wavelet  function  0(x)  and  the  magnitude  of  its  Fourier  spectrum  0(f) 
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filter,  possessing  both  symmetry  and  linear  phase.  The  wavelet  function  exhibits 
compact  support,  a  favorable  aspect  for  localizing  image  spectrum  content.  This 
condition  was  beneficial  for  reducing  aliasing  effects  brought  on  by  sampling  of  the 
wavelet  filter,  discussed  in  the  section  on  binary  detour  phase  holograms. 

3.1.2  Binary  Detour  Phase  Hologram.  Representing  the  complex- valued 
wavelet  filter  on  an  intensity-only  medium  (such  as  a  photographic  plate)  required 
holography.  A  binary,  detour  phase  hologram  represented  both  the  magnitude  and 
phase  of  the  complex  function  at  discrete  sample  intervals  (7).  The  first  consider¬ 
ation  when  constructing  the  hologram  was  the  extent  of  the  original  time  domain 
function.  According  to  sampling  theory,  the  wavelet  filter  must  be  sampled  at  small 
enough  intervals  in  order  to  prevent  aliasing.  The  sample  interval  must  not  exceed 
the  reciprocal  of  the  wavelet  function’s  extent.  In  this  manner,  the  compact  sup¬ 
port  of  the  wavelet  filter  helped  avoid  this  problem,  since  the  wavelet  function  had 
limited  extent.  The  detour  phase  hologram  was  divided  into  an  array  of  equal  di¬ 
mensional  cells.  Each  of  the  opaque  cells  corresponded  to  a  sample  point  in  the 
complex  function  plane.  Within  each  of  these  cells  was  placed  a  rectangular  aper¬ 
ture.  The  magnitude  and  phase  of  the  sample  point  defined  the  size  and  position  of 
the  aperture  within  the  cell.  Figure  6  shows  the  layout  for  one  cell. 

The  magnitude  of  the  sample  point  defined  the  height  of  the  rectangle:  hnm  = 
Anmdy.  The  rectangle  was  centered  vertically  within  the  cell  and  each  step  in  mag¬ 
nitude,  Anm,  equated  to  subsequent  increase  in  the  height  of  the  rectangle.  The 
magnitude  values  of  the  complex  function  were  scaled  (normalized)  so  that  the  max¬ 
imum  extent  of  the  function  corresponded  to  the  largest  extent  of  the  cell.  The 
width  of  each  rectangle,  w,  was  kept  constant.  The  value  of  the  width  determined 
the  efficiency  of  the  hologram.  A  wider  value  allowed  more  light  to  pass  through  the 
hologram,  which  was  an  important  factor  in  this  experiment.  However,  there  was 
a  trade-off  between  hologram  efficiency  and  distortion  due  to  phase  quantization. 
The  position  of  the  rectangular  aperture  was  determined  by  the  phase  of  the  sample 
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point,  cnm  =  Centering  the  aperture  equated  to  zero  phase,  shifting  to  the  right 
side  of  the  cell  was  +7r  phase,  and  shifting  the  rectangle  to  the  far  left  of  the  cell  was 
—7 r  phase.  The  number  of  quantization  levels  for  the  phase  was  determined  by  the 
ratio  of  the  cell  width  to  the  rectangle  width.  This  was  a  major  consideration  since 
the  primary  cause  of  distortion  of  the  original  image  by  this  holographic  technique 
lay  in  the  sampling  of  the  phase  (32).  The  higher  the  number  of  levels  the  lower 
the  distortion  caused  by  the  sampling  of  the  phase.  This  effect  was  best  illustrated 
by  the  relative  intensity  value  of  the  mth  order  false  image  to  the  zero  order  image: 

( mN  +  l)"2  (9) 

The  value,  TV,  in  equation  9  was  the  number  of  quantization  levels.  Increasing  N 
diminished  the  intensity  of  the  false  images,  and  limited  the  distortion.  The  holo¬ 
grams  for  this  research  were  created  using  a  laser  printer.  The  negative  image  was 
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transcribed  onto  a  transparency  to  be  used  to  created  a  photo-reduced  glass  slide. 
The  laser  printer  was  limited  to  a  400  dots-per-inch  resolution,  so  for  an  8  inch 
transparency,  3200  points  was  the  upper  limit  for  each  dimension.  The  holograms 
had  a  total  of  161  sample  points  in  each  dimension,  with  each  cell  made  up  of  16 
X  16  pixel  points.  The  161  sample  points  was  well  within  the  limit  set  by  sampling 
theory.  With  these  dimensions,  the  hologram  was  made  up  of  a  total  of  2576  points 
across  each  of  the  dimensions.  Because  the  actual  resolution  of  the  laser  printer 
fell  short  of  the  3200  point  limit,  this  value  seemed  appropriate.  The  width  of  each 
rectangle  was  chosen  to  be  8,  half  the  width  of  the  cell,  for  a  50  percent  duty  cycle. 
This  width  assured  that  the  hologram  would  pass  a  sufficient  amount  of  light  and 
still  allow  for  9  distinct  horizontal  positions  for  the  rectangle  within  each  of  the  cells. 
Nine  levels  of  phase  quantization  limited  the  distortion  according  to  equation  9.  The 
relative  intensity  for  the  first  order  false  image  was  one  one-hundredth  that  of  the 
original.  A  mathematical  analysis  of  the  array  of  cells  helped  to  understand  how 
the  size  and  position  of  the  these  rectangles  represented  magnitude  and  phase  of  a 
complex  function.  In  equation  10,  the  function  represented  the  hologram, 

which  was  simply  an  array  of  rectangular  aperatures.  For  this  analysis,  the  Fourier 
transform  of  the  array  of  rectangles  was  calculated  and  represented  by  the  function 
Ti(u,v).  This  function  was  rewritten  as  a  summation  over  the  array  dimensions  n 
and  m  as  shown,  which  is  the  discrete  Fourier  transform  of  the  original  function. 


T\(u,v) 


IT  (,(*.  y)e-i2’>"r*“‘>dxdy 

n  m 


(10) 


The  function  Gnm  was  the  Fourier  transform  of  the  rectangle  functions.  Each  Sine 
function  had  widths  defined  by  the  values  hnm  and  w,  the  height  and  width  of 
the  rectangle.  The  phase  factor  in  equation  11  was  due  to  the  shift  value  for  each 
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rectangle  within  the  cell,  c„m. 


Gnm(u,  v) 


sin  nwu  sin  irhnm v 
iru  irv 


(11) 


A  sample  point  equating  to  the  center  of  the  arrays  was  chosen  for  this  example. 
This  value  was  substituted  into  the  equation.  The  power  series  expansion  for  the 
Sine  function  was  evaluated  and  approximated  as  shown. 


Ue 


a. 


ve  =  0 


Gnm{u,  V )  =  sm  nwu hnme~i2wke'm/Jt  x  (1  +  nvh*m/6  -  j2ncnm(u  -  uc))  +  •  •  •  (12) 
ITU 

This  value  was  substituted  back  into  equation  10. 

Tl  (u,  „)  =  “E™  YTH  +  Qnm)hnme-*'ke"''d-  x  (13) 

™  n  m 

The  values  for  the  height  and  shift  were  substituted  back  into  the  equation  13: 
hnm  =  Anmdy  and  cnm  =  What  was  left  in  equation  14  was  the  discrete 

Fourier  transform  of  the  original  sample  point  of  magnitude  Anm  and  phase  <f>ntn. 
The  additional  factor  Qnm  was  due  to  phase  distortion.  This  analysis  showed  that 
the  array  of  rectangles  did  indeed  represent  the  original  complex  function. 

Tl(„,„)  =  yy  (i  +  Q„m)d,A„„e-‘k*~  x  (14) 

nU  n  m 

The  hologram  transparency  was  created  using  the  program  in  the  appendix.  The 
finished  product  was  a  series  of  opaque  rectangles  representing  the  negative  of  the 
hologram,  see  Figure  7.  This  transparency  was  used  to  create  a  photo-reduced  glass 
slide  for  use  on  an  optics  bench.  For  reducing  the  transparency,  the  negative  was 
taken  to  the  Cooperative  Electronics  Materials  Process  Lab  at  Wright  Patterson 
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AFB  for  reduction  onto  a  high  resolution  glass  plate.  The  Dekacon  optical  system 
used  for  reduction  of  the  CGH  is  shown  in  Figure  8. 

S.l.S  Two  Dimensional  Filters.  As  shown  in  Figure  1,  the  decomposition 
of  a  two  dimensional  signal  using  wavelet  analysis  resulted  in  combinations  of  the 
wavelet  and  basis  filters.  For  this  research  there  were  three  filters  required  to  capture 
the  vertical,  horizontal,  and  diagonal  details  of  the  image.  These  three  frequency 
channels  stemmed  from  the  combination  of  filtering  operations  across  the  rows  and 
columns  of  the  input  image.  In  this  manner,  combination  of  the  transforms  of  the 
wavelet  function  and  the  basis  function  resulted  in  horizontal  (Figure  9),  vertical 
(Figure  10),  and  diagonal  (Figure  11)  wavelet  filters.  Three  different  filters  were 
used  in  the  optical  experiment:  the  horizontal,  vertical,  and  diagonal  filter.  The 
filters  were  created  by  combination  of  the  basis  and  wavelet  function  depicted  in 
equations  7  and  8  respectively. 


Figure  7.  Transparency  of  Detour  Phase  Hologram 
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Figure  8.  Dekacon  Camera  System 


22 


Figure  9.  Horizontal  Wavelet  Filter:  <t>(x)ip(y) 
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Figure  10.  Vertical  Wavelet  Filter:  ip(x)<ff(y) 


Figure  11.  Diagonal  Wavelet  Filter: 
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3.2  Optical  Architecture 

In  order  to  first  characterize  the  optical  wavelet  transform,  a  MOSLM  was 
used  to  present  images  for  study.  The  optical  arrangement  was  developed  to  be  an 
improvement  on  the  MRA  algorithm.  The  MRA  algorithm  decomposed  the  input 
signal  into  its  wavelet  coefficients  at  dyadic  values  for  the  shift  and  scale  variable, 
necessary  for  swift  calculation  of  the  wavelet  transform.  Because  the  wavelet  trans¬ 
form  was  simply  a  correlation  between  the  input  image  and  the  wavelet  function, 
optics  implemented  the  wavelet  transform  in  real-time  with  continuous  distribution 
of  the  shift  variable.  This  was  an  advantage  over  the  MRA  algorithm.  Using  a  vari¬ 
ation  of  the  optical  correlator  allowed  continuous  distribution  of  the  scaling  function 
as  well. 

Continuous  distribution  of  the  scaling  function  was  realized  with  the  configu¬ 
ration  depicted  in  Figure  12.  This  was  an  improvement  over  the  optical  correlator 
in  Figure  3  because  the  negative  lens,  L2,  added  flexibility  over  the  scaling  of  the 
wavelet  filter.  In  the  configuration  in  Figure  12,  a  binarized  image  was  introduced  in 

Wavelet 

Input  Filter  Transform 

Plane  Plane  PJane 


the  input  plane  on  a  MOSLM.  Using  a  telescopic  lens  arrangement  to  perform  the 
Fourier  transforming  stages  in  the  optical  correlator  allowed  variation  of  the  scaling 
of  the  Fourier  transform.  This  arrangement  has  been  effective  in  developing  small 
scale  optical  correlators  and  can  now  be  used  for  scaling  the  wavelet  filter  (3).  The 
positive  lens,  LI,  was  placed  a  focal  length  away  from  the  input  plane.  The  Fourier 
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transform  of  the  image  appeared  at  the  back  focal  plane  of  lens  LI.  Lens  L2  imaged 
the  transformed  image  onto  the  filter  plane,  where  the  transform  was  a  virtual  object 
to  lens  L2.  As  long  as  the  transformed  image  fell  within  the  focal  length  of  L2,  a 
magnified  real  image  was  created  at  the  filter  plane.  The  placement  of  L2  relative 
to  the  back  focal  plane  of  LI  and  the  subsequent  placement  of  the  wavelet  filter  de¬ 
termined  the  scaling  of  the  transform  onto  the  wavelet  filter.  The  same  process  was 
duplicated  in  the  scaling  of  the  wavelet  transform  onto  the  camera  which  captured 
the  wavelet  transform  plane.  To  determine  the  minimum  scaling  of  the  wavelet  fil¬ 
ter,  the  system  was  analyzed  according  to  the  dimensions  of  the  MOSLM  and  the 
wavelet  filter. 

The  minimum  focal  length  for  the  transforming  lens  in  a  matched  spatial  filter 
is  determined  by  the  spacing  of  the  input  image  pixels  (di),  the  spacing  of  the  filter 
plane  pixels  (d2),  the  dimension  of  the  filter  (in  pixels)  (jV2),  and  the  wavelength  of 
the  laser  source  (A)  (11). 


F  =  (15) 

A 

With  the  MOSLM  pixel  spacing  of  76  /xm,  the  filter  plane  pixel  spacings  of 
63  /im,  the  filter  plane  pixel  dimension  of  161  element  array,  and  the  wavelength 
of  the  laser  at  632.8  nm;  the  minimum  focal  length  was  1.22  meters.  This  focal 
length  was  the  to  =  0  or  original  wavelet  for  this  experiment.  Repositioning  of 
lens  L2  and  the  filter  plane  allowed  a  continuous  variation  of  the  scaling  variable  by 
expanding  the  Fourier  transform  onto  the  wavelet  filter  plane.  In  this  arrangement, 
only  one  hologram  for  each  of  the  wavelet  filters  was  needed  to  produce  multiple 
scale  resolutions. 
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Figure  13.  Improved  Optical  Wavelet  Fingerprint  Identification  Arrangement 
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3.3  Optical  Wavelet  Fingerprint  Identification 

After  the  wavelet  filter  was  characterized,  it  was  used  in  an  optical  fingerprint 
identifier.  The  optics  bench  arrangement  is  shown  in  Figure  13.  In  this  arrangement 
a  fingerprint  image  was  scanned  into  the  system  using  a  90  degrees  prism.  By  placing 
the  person’s  finger  onto  the  prism,  the  areas  where  the  fingerprint  ridges  came  in 
direct  contact  with  the  surface  of  the  prism  were  areas  where  the  laser  beam  was 
absorbed  and  not  reflected  into  the  system.  Thus,  the  fingerprint  image  was  read 
into  the  system  using  frustrated  total  internal  reflectance. 

The  optical  fingerprint  identifier  used  a  matched  spatial  filter  (MSF)  for  dis¬ 
criminating  between  fingerprint  images.  In  order  to  create  the  MSF,  the  second 
hologram  was  created  using  Vander  Lugt  filtering. 


P 


Figure  14.  Vander  Lugt  Filter  (18). 

3.3.1  Vander  Lugt  Filter.  The  MSF  was  created  through  use  of  Vander- 
Lugt’s  method  of  interferometrically  recorded  frequency-plane  masks  (22).  The 
optical  arrangement  is  shown  in  Figure  14.  The  Fourier  spectrum  of  the  wavelet 
transformed  fingerprint  image  was  realized  in  the  back  focal  plane  of  the  lens,  L3  in 
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Figure  13.  A  reference  beam  was  split  off  from  the  laser  beam  and  interfered  with  the 
spectral  content  of  the  wavelet  transformed  fingerprint  image.  In  this  manner  the 
phase  of  the  complex  field  was  saved  through  the  interference  of  the  plane  wave  and 
the  complex  field.  The  resulting  interference  field  was  captured  using  a  thermoplastic 
plate.  The  HC-300  Holographic  recording  camera  was  used  to  store  the  interference 
pattern  of  the  Vander  Lugt  filter.  The  HC-300  developed  the  MSF  electrically  in 
about  one  minute.  The  thermoplastic  plate  was  developed  using  a  four  step  process 
shown  in  Figure  15.  After  the  plate  was  developed,  it  was  used  as  a  MSF  against 
which  other  wavelet  transformed  images  could  be  compared.  The  fingerprint  image 
was  correlated  with  the  wavelet  filter,  then  the  wavelet  transform  of  the  fingerprint 
was  correlated  with  the  MSF.  The  energy  in  the  correlation  peak  captured  in  the 
diffracted  beam  off  the  hologram  of  the  Vander  Lugt  filter  was  used  to  discriminate 
between  fingerprints. 
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THERMOPLASTIC 
PHOTOCONDUCTOR 


a.  CHARGING 


k  uniform  chir|t  is  deposited  on  thi  thermoplastic  tad  photocoaductor. 


Exposure  to  laser  flings  pattern  redistributes  the  charges  through 
the  photoconductor. 


c.  SECOND 
CHARGING 


♦++  +++ 
+++++-M-+++ 


Recharge  increases  the  electric  field  across  the  exposed  area  of 
the  thermoplastic. 


d.  DEVELOPMENT 


Hooting  of  the  thermoplastic  causes  permanent  deformation.  The 
irregularities  will  defract  light  to  recreate  the  original  image. 


Figure  15.  HC-500  Thermoplastic  Recording  Process  (24). 
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IV.  Results 


This  chapter  presents  and  discusses  the  results  based  on  the  methodology  de¬ 
scribed  in  Chapter  III.  The  two  topics  of  this  research,  application  of  an  optical 
biorthogonal  wavelet  transform  and  development  of  an  improved  optical  wavelet 
fingerprint  identifier,  are  examined.  Section  4.1  compares  the  optical  biorthogonal 
wavelet  transform  to  a  digital  wavelet  transform  of  an  image  of  the  letter  “E” .  Sec¬ 
tion  4.2  shows  results  from  an  optical  biorthogonal  wavelet  transform  of  an  image 
of  “Lenna”  and  a  fingerprint  at  three  resolutions.  The  results  from  an  improved 
optical  wavelet  fingerprint  identifier  are  shown  in  Section  4.4.  Finally,  the  rotational 
invariance  of  the  improved  wavelet  optical  fingerprint  identifier  are  compared  to  the 
standard  Vander  Lugt  matched  filter  method. 

4-1  Characterization  of  Optical  Biorthogonal  Wavelet  Transform 

An  image  of  the  letter  “E”  was  used  in  characterizing  the  optical  process.  The 
letter  “E”  was  chosen  because  it  contained  easily  recognized  horizontal,  vertical,  and 
corner  (diagonal)  features.  Three  “types”  of  features  were  required  to  accurately 
characterize  the  set  of  three  wavelet  filters  shown  in  Figures  9-11. 

4-1.1  Input  Image  of  Letter  “E”.  The  image  of  the  letter  “E”  is  shown 
in  Figure  16.  The  original  image  was  displayed  on  a  MOSLM  which  had  an  array 
of  128  X  128  pixels.  The  image  shown  in  Figure  16  was  captured  at  the  output 
of  the  optical  correlator  shown  in  figure  3  with  the  wavelet  filter  removed  from  the 
setup.  Perfect  alignment  of  the  four  lenses  was  assured  by  examining  the  output 
image  shown  here. 

4-1.2  Horizontal  Filter  Results.  The  results  from  the  horizontal  filter  are 
shown  in  Figures  17  through  20.  The  optical  transform  for  the  finest  resolution 
horizontal  wavelet  filter  is  shown  in  Figure  17.  This  corresponds  to  level  m  =  0  in 
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Figure  16.  Capital  Letter  “Eu 


equations  3  and  4.  The  digitally  produced  transform  for  the  same  resolution  wavelet 
filter  is  shown  in  figure  18.  As  shown  in  the  figures,  the  horizontal  structure  of  the 
image  was  emphasized  by  the  horizontal  filter.  The  horizontal  details  have  been 
isolated  from  the  rest  of  the  image.  The  results  from  the  optical  process  compared 
favorably  with  those  created  digitally. 

The  optical  and  digital  results  for  the  first  dyadic  scaling  of  the  wavelet  are 
displayed  in  Figures  19  and  20  respectively.  As  the  wavelet  function  was  doubled  in 
size,  its  filter  was  diminished  the  same  amount.  The  result  was  a  filter  which  passed 
a  lower  band  of  spatial  frequencies.  Less  horizontal  structure  was  observed  in  the 
transformed  image  as  expected.  Again,  the  optical  results  were  very  similar  to  the 
digitally  produced  transform. 
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Figure  17.  Optical  Wavelet  Transform  for  Horizontal  Filter  on  Letter  E,  m  =  0 
Resolution 


Figure  18.  Digital  Wavelet  Transform  for  Horizontal  Filter  on  Letter  E,  m  =  0 
Resolution 
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Figure  20.  Digital  Wavelet  Transform  for  Horizontal  Filter  on  Letter  E,  m  =  1 
Resolution 
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4- 1-3  Vertical  Filter  Results.  The  optical  transform  for  the  finest  resolution 
vertical  wavelet  filter  is  shown  in  Figure  17,  and  the  digitally  produced  transform  for 
the  same  wavelet  filter  is  shown  in  Figure  18.  The  optical  and  digital  results  for  the 
first  dyadic  scaling  of  the  wavelet  are  displayed  in  Figures  23  and  24  respectively. 
The  results  from  the  optical  process  again  compared  favorably  with  those  created 
digitally.  In  this  case,  the  vertical  details  were  isolated  by  the  wavelet  transform. 

4-1-4  Diagonal  Filter  Results.  The  optical  transform  for  the  finest  res¬ 
olution  diagonal  wavelet  filter  is  shown  in  Figure  25,  and  the  digitally  produced 
transform  for  the  same  wavelet  filter  is  shown  in  figure  26.  The  optical  and  digital 
results  for  the  first  dyadic  scaling  of  the  wavelet  are  displayed  on  Figures  27  and 
28  respectively.  For  the  diagonal  filter,  the  results  from  the  optical  process  did  not 
compare  well  with  those  of  the  digital  ones.  The  images  were  faint  and  distorted  in 
comparison  to  the  other  filters.  This  was  due  to  the  input  image.  Reviewing  the 
image  produced  by  the  MOSLM  in  Figure  16,  the  comers  of  the  letter  were  less  than 
optimum.  The  MOSLM  was  old  and  missing  some  of  the  pixels  in  the  array,  making 
the  image  distorted.  The  lack  of  distinction  in  the  comers  reduced  the  amount  of 
energy  in  those  spatial  frequencies  which  were  passed  by  the  diagonal  wavelet  filter. 
The  limited  intensity  of  the  optical  wavelet  transform  was  due  to  this  distortion  in 
the  input  image. 
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Figure  21.  Optical  Wavelet  Transform  for  Vertical  Filter  on  Letter  E, 
Resolution 


Figure  22.  Digital  Wavelet  Transform  for  Vertical  Filter  on  Letter  E, 
Resolution 


Figure  23.  Optical  Wavelet  Transform  for  Vertical  Filter  on  Letter  E, 
Resolution 


Figure  24.  Digital  Wavelet  Transform  for  Vertical  Filter  on  Letter  E, 
Resolution 


Figure  25.  Optical  Wavelet  Transform  for  Diagonal  Filter  on  Letter  E, 
Resolution 


Figure  26.  Digital  Wavelet  Transform  for  Diagonal  Filter  on  Letter  E, 
Resolution 


Figure  27.  Optical  Wavelet  Transform  for  Diagonal  Filter  on  Letter  E,  m  =  1 
Resolution 
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Figure  28.  Digital  Wavelet  Transform  for  Diagonal  Filter  on  Letter  E,  m  =  1 
Resolution 
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4.2  Optical  Biorthogonal  Wavelet  Transform  of  henna 

The  image  of  Lenna  was  used  in  order  to  compare  to  results  obtained  in  past 
research  efforts  at  the  Air  Force  Institute  of  Technology.  The  image  of  Lenna  ap¬ 
pears  in  this  research  because  it  contains  significant  complexity.  The  image  contains 
detail  at  various  resolution,  optimal  for  exploitation  using  multiresolutional  analysis. 
Past  research  efforts  explored  correlation  with  a  Vander  Lugt  Matched  Spatial  Filter 
(MSF)  of  a  Harr  wavelet  (25,  8)  as  well  as  correlation  with  a  phase-only  spatial  filter 
of  the  Harr  wavelet  created  using  a  computer  generated  hologram  (2). 

The  image  of  Lenna  was  binarized  and  placed  on  the  MOSLM.  As  before  with 
the  letter  “E  “,  the  image  of  Lenna  was  captured  at  the  output  of  the  optical  correlator 
shown  in  figure  3  with  the  wavelet  filter  removed  from  the  setup.  The  image  is  shown 
in  Figure  29.  The  image  of  Lenna  was  useful  for  exploring  the  wavelet  transform 
because  of  the  various  levels  of  detail  which  exist  at  different  scales  in  the  image. 
Each  of  the  three  wavelet  filters  depicted  in  Figures  9-11  was  applied  to  the  image 


Figure  29.  Lenna 


of  Lenna  with  good  results. 
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Each  of  the  three  wavelet  filters  was  applied  to  the  image.  Results  showed 
how  the  various  image  details  were  enhanced  by  each  of  the  filters.  In  turn,  three 
resolutional  scalings  of  each  wavelet  were  applied  to  the  image  of  Lenna,  through 
adjustment  of  the  optical  correlator  configuration.  This  resulted  in  nine  different 
wavelet  transforms  of  the  image  of  Lenna. 

The  three  scalings  of  the  horizontal  wavelet  are  shown  in  Figure  30.  The  finest 
resolution  horizontal  wavelet  represented  the  higher  passband  of  spatial  frequencies. 
The  edges  of  the  image  represented  the  highest  frequency  components.  The  eyebrows, 
the  mouth,  and  the  brim  of  the  hat  were  enhanced  at  this  resolution.  As  the  scaling 
factor  was  doubled,  the  effect  was  to  allow  some  of  the  lower  spatial  frequencies  to 
pass.  More  details  from  the  hat  and  face  were  visible  now.  Other  objects  in  the 
background  became  apparent.  As  the  scaling  was  doubled  again,  the  image  details 
began  to  take  on  more  form.  The  lines  have  been  replaced  with  the  low  frequency 
components  making  up  the  shapes  within  the  image.  The  three  step  transition 
between  linear  details  and  form  was  obvious  with  the  vertical  wavelet  as  well.  In 
Figure  31  the  vertical  details;  the  edges  of  the  nose,  hair,  and  feathers,  were  the 
target.  Again,  as  the  scale  of  the  wavelet  was  increased,  more  of  the  form  became 
more  prominent  as  the  edges  were  left  out.  Finally  with  the  diagonal  wavelet,  in 
Figure  32  the  process  was  repeated.  This  time  the  shoulder,  brim  and  crown  of  the 
hat,  were  enhanced  in  the  image.  Each  of  these  nine  images  represented  a  unique 
resolutional  detail  for  the  image  of  Lenna.  It  was  obvious  from  these  results  that 
the  optical  wavelet  could  be  applied  to  the  fingerprint  to  extract  similar  features. 
Features  which  could  be  used  in  the  identification  process. 

4.3  Optical  Biorthogonal  Wavelet  Transform  of  a  Fingerprint 

The  same  nine  wavelet  transforms  were  applied  to  a  fingerprint  image.  This 
process  revealed  several  levels  of  detail  within  the  fingerprint,  detail  which  could 
be  exploited  for  recognition  algorithms.  The  input  image  of  the  fingerprint  was 
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Figure  30.  Horizontal  Wavelet  Transform  Results  for  Lenna:  (a)  m  =  0  (b)  m  =  1 
(c)  m  =  2 
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(C) 

Figure  31.  Vertical  Wavelet  Transform  Results  for  Lennar  (a)  m  —  0  (b)  m  =  1  (c) 
m  =  2 
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Figure  32.  Diagonal  Wavelet  TVansform  Results  for  Lenna:  (a)  m  =  0  (b)  m  =  1 
(c)  m  =  2 
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binarized  and  placed  on  the  MOSLM.  The  image  of  the  fingerprint  as  captured 
through  the  CCD  camera  is  shown  on  Figure  33.  The  fingerprint  image  consisted 


Figure  33.  Fingerprint 


mainly  of  ridge  lines.  Since  it  was  the  structure  of  these  lines  which  uniquely  specified 
an  individual’s  fingerprint,  the  wavelet  filter  was  used  to  enhance  the  structure  of 
these  lines  at  various  spatial  frequencies.  The  fingerprint  image  was  almost  radial 
in  structure,  containing  lines  which  had  spatial  frequencies  in  all  three  locations: 
horizontal,  vertical,  and  diagonal.  As  before  with  the  image  of  Lenna,  each  of  the 
nine  wavelet  transformed  images  was  unique.  Of  special  interest  was  the  results  with 
the  finest  resolution  diagonal  wavelet  transform  in  Figure  36(a).  In  this  case  there 
was  plenty  of  energy  in  the  diagonal  frequencies,  avoiding  the  problem  experienced 
earlier  with  the  letter  “E”  (see  Figure  25). 

To  enhance  the  optical  fingerprint  identification  process,  the  digital  wavelet 
transform  was  applied  to  the  fingerprint.  The  optical  correlator  identification  pro¬ 
cess  was  limited  by  its  rotational  variance.  To  remedy  this,  the  wavelet  transform 
was  used  to  enhance  those  features  of  the  fingerprint  that  were  rotational  symmet- 
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rical.  Since  the  fingerprint  was  mainly  radial  in  structure,  the  wavelet  transform 
could  be  used  to  emphasize  that  radial  structure  for  identification.  The  finest  reso¬ 
lution  diagonal  wavelet  transform  was  interesting  because  of  two  reasons.  First,  it 
contained  a  significant  amount  of  energy.  Energy  which  was  required  for  creating 
the  thermoplastic  hologram  in  the  identification  stage.  Second  and  more  important, 
the  diagonal  wavelet  transform  enhanced  the  radial  line  structure  in  a  way  which 
approached  rotational  symmetry.  This  transform  was  chosen  as  the  best  candidate 
for  the  improved  optical  wavelet  fingerprint  identifier. 
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(C) 


Figure  34.  Horizontal  Wavelet  Transform  Results  for  Fingerprint:  (a)  m  =  0  (b) 
m  =  1  (c)  m  =  2 
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4  4  Improved  Optical  Wavelet  Fingerprint  Identifier 

A  set  of  subjects  was  used  for  identification,  the  same  set  used  in  parallel  per¬ 
sonal  identification  research.  These  other  research  efforts  involved  speech,  handwrit¬ 
ing,  and  face  recognition  (26,  28,  19).  The  subjects  placed  their  thumbprint  on  the 
prism  in  the  optical  arrangement  portrayed  in  Figure  13.  The  computer  generated 
hologram  of  the  diagonal  wavelet  filter  was  used  to  enhance  the  diagonal  features 
of  the  fingerprint  by  computing  the  wavelet  transform  of  the  print.  This  image  was 
than  correlated  with  a  Vander  Lugt  MSF  of  a  wavelet  transformed  print.  The  corre¬ 
lation  peak  was  then  analyzed  to  determine  if  this  arrangement  could  discriminate 
between  wavelet  transformed  fingerprints. 

Results  showed  that  the  wavelet  transformed  fingerprint  did  indeed  contain 
enough  uniqueness  to  discriminate  between  fingerprints.  The  correlation  peaks  cap¬ 
tured  with  the  CCD  camera  are  shown  in  Figure  37.  A  thermoplastic  hologram 
was  created  using  the  first  subject’s  print.  The  auto-correlation  peak  was  captured 
with  the  print  left  on  the  prism  after  creating  the  thermoplastic  hologram.  Subject 
one  then  removed  the  thumb  and  subsequently  returned  the  thumb  to  the  prism, 
after  it  had  been  wiped  clean.  The  resulting  correlation  peak  was  then  captured 
and  recorded  as  “subject-one  correlation  “.  Four  other  subjects  were  than  used  to 
match  against  subject  one.  The  result  was  that  the  correlation  peak  was  much  more 
obvious  for  subject  one.  This  proved  that  the  wavelet  transformed  fingerprint,  with 
its  reduced  feature  set  could  still  discriminate  between  fingerprints. 
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Auto-correlation  Subject  1  Cross-correlation 

With  Subject  1  MSF 


Subject  2  Cross-correlation  Subject  3  Cross-correlation 

With  Subject  1  MSF  With  Subject  1  MSF 


Subject  4  Cross-correlation  Subject  5  Cross-correlation 

With  Subject  1  MSF  With  Subject  1  MSF 


Figure  37.  Improved  Optical  Wavelet  Fingerprint  Identifier  Results:  Discrimina 
tion  Against  Subject  1  MSF 
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4-5  Test  for  Rotational  Invariance 

To  determine  if  the  wavelet  transformed  fingerprint  added  any  rotational  in¬ 
variance  to  the  Vander  Lugt  MSF,  a  comparison  was  made  between  the  standard 
4Vander  Lugt  MSF  correlator  and  the  wavelet  process.  A  glass  slide  of  a  fingerprint 
was  used  as  the  input  image.  The  glass  slide  was  placed  on  a  rotating  stand  which 
could  measure  the  angular  rotation  in  degrees.  A  MSF  was  created  for  both  the 
normal  and  wavelet  transformed  non- rotated  fingerprint.  The  auto-correlation  peak 
was  recorded.  Then  the  glass  slide  with  the  print  on  it  was  rotated  so  that  the  effect 
of  rotation  could  be  observed  and  recorded.  The  results  are  shown  in  Figure  38, 
where  the  correlation  peaks  were  captured  for  the  normal  and  the  wavelet  trans¬ 
formed  fingerprint.  The  correlation  peaks  for  the  offset  results  are  plotted  in  Figures 
40  through  45.  The  wavelet  fingerprint  displayed  a  greater  amount  of  invariance  to 
the  rotation  out  to  5  degrees  off  center.  The  energy  in  the  peak  remained  strong 
in  direct  contrast  with  the  peak  for  the  normal  fingerprint.  The  wavelet  correlation 
peak  was  strong  enough  in  the  5  degrees  offset  so  that  discrimination  was  possible 
between  the  rotated  fingerprint  and  a  different  on-axis  print.  To  illustrate  this,  the  5 
degree  offset  peak  was  compared  to  the  results  in  Figure  37  for  subject  2.  As  shown 
in  Figure  39,  the  offset  fingerprint  was  still  distinctive  from  another  subject.  Even 
when  the  fingerprint  was  offset  by  as  much  as  five  degrees,  the  optical  identifier  was 
able  to  determine  the  identity.  This  adds  flexibility  to  the  optical  identification  pro¬ 
cess,  allowing  the  subject  some  latitude  when  placing  the  print  in  the  arrangement 
for  identification. 

Comparison  of  the  signal-to-noise  ratio  (SNR)  for  the  correlation  peaks  in 
Figure  38  is  shown  in  Table  2.  The  value  for  the  SNR  was  determined  by  comparing 
the  average  value  of  the  signal  to  that  of  the  noise.  To  discriminate  between  signal 
and  noise,  the  largest  pixel  value  was  determined.  The  signal  was  considered  to  be 
those  pixel  values  within  3  dB  of  the  peak.  Everything  else  was  considered  to  be 
noise. 
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Normal  Autocorrelation  Wavelet  Autocorrelation 


3  Degree  Offset 


Normal  Autocorrelation 


Wavelet  Autocorrelation 


Normal  Autocorrelation  Wavelet  Autocorrelation 

Figure  38.  Auto-Correlation  Peak  vs  Offset  Rotation  for  Normal  and  Wavelet 
Transformed  Fingerprint  Image  at  Various  Rotations 
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Offset 

SNR  Wavelet  Correlator 

SNR  Normal  Correlator 

0° 

51.34 

49.72 

K9 

38.71 

22.98 

R9 

40.5 

20.74 

Table  2.  Signal-to-Noise  Ratios  for  Various  Offsets 


The  values  shown  in  Table  2  were  generated  by  comparing  the  average  value 
of  the  signal  peak  to  the  standard  deviation  of  the  noise.  As  shown  in  equation  16, 
the  average  value  of  the  noise  was  subtracted  from  the  peak  in  order  to  reduce  bias 
(20). 


SNR  = 


fi(Peak)  —  ji(Noise) 
ff(Noise) 


(16) 


(a)  (b) 

Figure  39.  Comparison  of  (a)  Five  Degree  Offset  Wavelet  Auto-correlation  and  \h) 
Non-Rotated  Subject  Two  Wavelet  Cross-correlation  with  Subject  One 
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Figure  42.  MSF  Correlation  Peak  Plot  for  Normal  3  Degree  Offset  Fingerprint 


Figure  43.  MSF  Correlation  Peak  Plot  for  Wavelet  Transformed  3  Degree  Offset 
Fingerprint 


Figure  44.  MSF  Correlation  Peak  Plot  for  Normal  5  Degree  Offset  Fingerprint 


Figure  45.  MSF  Correlation  Peak  Plot  for  Wavelet  Transformed  5  Degree  Offset 
Fingerprint 
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4-6  Summary 

This  chapter  showed  the  results  of  the  biorthogonal  wavelet  transform  on  three 
input  images.  The  results  were  compared  to  digital  simulations  and  past  research 
efforts.  The  wavelet  transform  was  used  to  improve  upon  the  optical  fingerprint  iden¬ 
tifier.  The  results  were  favorable.  SNR  ratios  for  the  wavelet  transformed  fingerprint 
showed  a  greater  tolerance  to  rotational  offsets  than  did  the  normal  fingerprint.  The 
next  chapter  draws  conclusions  and  provides  suggestions  for  future  research. 
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V.  Conclusions 


This  research  implemented  the  optical  wavelet  transform  for  a  symmetric  biorthog- 
onal  wavelet  function.  The  optical  wavelet  transform  demonstrated  distinct  advan¬ 
tage  over  the  computer  based  MRA  algorithm.  The  optical  process  computed  the 
wavelet  transform  in  real-time  with  continuous  shift  variable.  The  scale  variable  of 
the  wavelet  transform  was  varied  with  continuous  distribution  through  manipulation 
of  optical  arrangement.  This  was  an  advantage  over  the  computer  based  technique 
which  uses  discrete,  dyadic  values  for  scale  and  shift  variables. 

A  detour  phase  hologram  of  the  two  dimensional  wavelet  filters  was  constructed 
and  operational.  The  optical  wavelet  transformed  images  compared  well  with  digital 
methods.  This  optical  method  also  showed  improvement  over  past  research  in  this 
area  here  at  AFIT.  The  choice  of  the  holographic  technique  and  type  of  wavelet 
showed  great  promise. 

The  optical  wavelet  transform  of  a  fingerprint  added  to  the  identification  ro- 
cess.  An  optically  transformed  fingerprint  image,  a  reduced  feature  set  of  the  origi¬ 
nal,  still  demonstrated  enough  unique  qualities  for  discrimination  with  an  MSF.  The 
wavelet  transformed  image  also  demonstrated  an  added  rotational  invariance  when 
compared  to  a  normal  fingerprint  image.  This  successful  application  of  the  optical 
wavelet  transform  gives  hope  to  further  study. 

5.1  Future  Applications 

Using  the  optical  wavelet  transform  for  enhancing  fingerprint  identification 
would  assist  in  the  development  of  a  field  unit  for  scanning  fingerprints.  This  ap¬ 
plication  would  benefit  law  enforcement  agencies,  allowing  rapid  storage,  identifi¬ 
cation,  and  retrieval  of  fingerprint  images,  with  the  rotational  invariance  allowing 
some  degree  of  flexibility  when  attempting  to  input  a  print.  The  enhanced  fingerprint 
identification  would  also  improve  security  access  systems  in  remote  locations. 
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The  wavelet  transform  also  shows  promise  for  image  segmentation.  An  optical 
process  could  be  used  in  automatic  pattern  recognition  for  application  in  a  smart 
weapon  system.  A  system  could  be  developed  which  could  effectively  scan  for  several 
resolutions,  using  those  resolutions  to  discriminate  between  objects  in  a  scene. 

5.2  Recommendations 

For  continued  research  in  optical  wavelet  transforms,  the  following  areas  should 
be  further  investigated: 

1.  A  comparison  of  different  wavelet  functions  for  transforming  fingerprint 
images  should  be  made  to  understand  how  the  choice  of  wavelet  will  affect  the 
identification  process. 

2.  Investigate  the  use  of  the  optical  wavelet  transform  in  image  segmentation 
for  pattern  recognition. 

3.  Demonstrate  an  optical  process  for  decomposition  and  perfect  reconstruction 
of  an  image. 

4.  A  comparison  of  other  magnitude  and  phase  holographic  techniques  for  cre¬ 
ating  the  wavelet  filter  would  determine  the  optimal  method  for  limiting  distortion. 

5.  And  finally,  access  to  a  larger  printer  with  similar  resolution  would  allow  a 
larger  amount  of  phase  quantization  levels,  limiting  the  major  cause  of  distortion  in 
the  holographic  process. 
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Appendix  A.  C  Program  used  to  Create  Detour  Phase  Computer 

Generated  Hologram 

•include  <stdio.h> 

•include  <math.h> 

•include  "macros. h" 

float  *vector(); 
float  **matrix(); 

main(argc.argv) 
int  argc ; 
char  **argv; 

{int  x,y,i,j,  temp,  bxsz,  tempi,  temp2,  bxht,  bxvd; 
float  *rdata; 
float  **mac; 

double  pi,  hmag,  **magdata,  *mag,  **phasedata,  *phase; 
int  number. cells, numbe r _po int s ; 

FILE  *magf ile,  *phasefile,  *output_f ile,  *magchk,  *phschk; 

number_points=2576 ; 
number_cells=161 ; 
pi=3. 1516152654; 

/*  Allocate  memory  for  the  input  files  for  magnitude  and  phase  of  vavelet  filter  */ 

magdata=dmatrix(l .number.cells, 1 .number.cells) ; 
mag=dvector(l .number.cells) ; 

phasedata=dmatrix(l .number.cells , 1 ,number_cells) ; 
phase=dvector(l ,number_cells) ; 

rdata=vector(l ,number_points) ; 

mac=matrix(l .number .points , 1 .number. points) ; 

/*  Output  file  for  CGH  data  */ 

output.file  *  fopenC'cgh"  ,"w"); 


n&gf i le*f open ( "magcghps i " ,  "r") ; 
phasef ile*f openO'phasecghpsi" ,  "r") ; 

■agchk=f openC'nagchk" , "v") ; 
phschk=f  openO'phschk" , "w") ; 

loopli (number  cells) 

{ 

fscanf (magf ile,  "Xle",  taagCi]); 
fscanf (phasef ile ,  "Xle",  kphaseCi]); 

> 

loopli (number.cells) 

{ 

loopl j (number.cells) 

{ 

magdata[i]  [j]=mag[i] *mag [j]  ; 
phasedataCi]  [j]=phase[i]+phase[j]  ; 

fprintf (magchk,  "Xd  \n",  2*Cint)(magdata[i] Cj]*42.44/2)) ; 

if  (phasedataCi] [j]>pi) 

{phasedataCi] [j]=phasedata[i] [j]-(2*pi) ;> 
if  (phasedata[i] [j]<-pi) 

{phasedataCi]  Cj]=phasedataCi]  Cj]+(2*pi)  ;> 

fprintf (phschk,  "Xd  \n",  (int) (floor (phasedataCi] Cj]*4/pi))) ; 

> 

> 

for  (y=l ; y<=number_cells ; y++) 

{ 

for  (x=l;x<=number_cells;x++) 

{ 


{templ=2*(int) (magdataCx] Cy]*42.44/2) ;} 
{bxht=(int)((16  -  templ)/2);> 
if  (phasedata  Cx] Cy] <0) 

{temp2=(int) (floor (phasedata Cx] Cy]*4/pi)) ;> 
else 
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{temp2*(int)(ceil(phasedata[x] Cy]*4/pi)) ;> 
{bxvd*(int)((6  -  temp2)/2) ;} 


for  (i  *  16*y  -  15;i<*16*y;i++) 

{for  (j*  16*x  -  IS; j<*16*x; j-*-0 

{ 

if  (i>«16*y-15+bxht  k  i<16*y-15+bxht+templ  k  i>*16*x-ll+temp2  k  i<16*x-3+temp2 
{mac [j] [i]*1.0;> 

else 

{mac[j]  [i]*0.0;> 

>  /*  end  of  j  loop  */ 

>  /*  end  of  i  loop  */ 


>  /*  end  x  loop  */ 

>  /*  end  y  loop  */ 


for  (i  =  1908; i<*1956;i++) 

{for  (j*  1280 ;j <*1296 ;j++) 
{mac [i]  [j]=l  .0; 

> 

> 

for  (j  *  1908; j <=1956; j++) 

{for  (i=  1280 ; i<=1296 ; i++) 
{mac [i] [j]*1.0; 

> 

> 


printf ("Values  computed  saving  the  data\n"); 
ff lush(stdout) ; 

for  (y*l ; y<=number_points ; y++) 

{ 

for  (x=l ;x<=number_points;x++) 

{rdata [x]* (float) mac [x] [y] ;> 

f write (rdata,sizeof (float) .number^points, output _f ile) ;> 
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free_vector(rdata,  1,  nuaber .points) ; 
free.aatrixCaac,  1,  nuaber .points,  1,  nuaber.points) ; 
free.dvectorCaag,  1,  nuaber.cells); 
free.dvector (phase,  1,  nuaber.cells); 
free.daatrixCaagdata,  1,  nuaber.cells,  1,  nuaber.cells); 
free.daatrixCphasedata,  1,  nuaber.cells,  1,  nuaber.cells); 
f close (output.f ile) ; 
fclose(magf ile) ; 
f close (phasef ile) ; 
f close (magchk); 
f close (phschk) ; 

>  /*  end  main  */ 
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